Differential abundance analysis using ANCOMBC was carried out on unrarefied data due to the recommendations by the package authors and discussions online. For significant taxa associated with different site-treatment groups, the data was non-aggregated to ensure there were enough samples for ANCOM analysis.
No taxa at phyla or genera level were identified as differentially abundant when comparing control and water stress samples overall.
We also tested for phyla, families or genera which were differentially abundant between control and drought on individual sites. Again, ANCOMBC did not identify any significant taxa. This was surprising given PERMANOVA tests determined that treatment had a significant effect on community composition for unweighted UniFrac.
A reason for this could be the small sample sizes. Bioinformatics tools such as ANCOMBC are less effective for low sample counts, while small sample sizes can increase the risk of Type II errors. These issues may have hampered our ability to find significant differences between treatments
See further down for results which were not significant
However, when we compared sites within treatments, we had some significant results. Multiple phyla and genera passed both ANCOMBC’s differential abundance and sensitivity tests.
For both control and drought, several phyla appeared as differentially abundant between sites in both treatments. For example, Gemmatimonadota, Chloroflexota and Eisenbacteria.
Some phyla were not differentially abundant in control samples but showed strong differences between sites in the drought samples, indicating a drought response.
Below are the results of the global test for controls. Site comparisons (in the tabs) have been split into separate tables containing all comparisons for that particular site.
For our analysis, we only considered taxa which had passed both the differential abundance and sensitivity tests. Taxa which passed one test but not the other were disregarded. However, the results from all taxa are shown in the tables below.
For the control samples, differences in abundance under ambient conditions indicated significant phyla were responding to strongly to site characteristics.
Phyla showing the strongest differences between sites in the controls (and which passed both tests) were Nitrospirota A and Gemmatimonadota.
Nitrospirota A contain nitrifying taxa, and is a phyla which is generally unresponsive to drought, but does respond to site factors such as land management intensity and soil nitrogen availability.
For the drought samples, Bacillota A-368345, which was not significantly different between sites’ control samples, showed a very strong change in abundance between sites. Bacillota A-368345, a strain in the phyla formerly known as Firmicutes, are drought resistant. Their ability to withstand or thrive in drought conditions is linked to their capacity to form endospores, which remain dormant in the soil and confer a competitive advantage in droughted-soils
Site pairwise comparisons is where we could see several changes happening. Major changes included:
Due to time contraints, we were not able to carry out a detailed analysis of changes at a genera level, but there were some significant results for both control and drought samples. Full tables are below
As with bacteria, no differentially abundant phyla associated with treatment, either overall or on individual sites, were identified by ANCOMBC analysis.
Again, this may be due to small sample sizes. The extra challenge we had was missing samples for some sites due to quality control failures. ANCOMBC needs at least three samples per factor and both Windmill Farm and Castle Field West had only two samples in some treatments. This meant, treatment comparisons on these sites could not be carried out.
When comparing sites within treatments, we had some significant results although not as many as for bacteria. Additionally, the p-values in the global tests were not as low as for bacteria, indicating a weaker response, particularly to drought conditions.
For control samples, only four phyla were differentially abundant in global tests:
For drought samples only one phylum was differentially abundant between sites: Zoopagomycota
Time constraints meant we were unable to complete a full analysis of how abundances change for specific genera. The results tables are below.
Below are the outputs from significance tests for treatment, where there were no significant results. These are included for academic transparency and in case they are of interest.